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.^ . Abstract In this survey article we present difference field algorithms for symbolic 

summation. Special emphasize is put on new aspects in how the summation prob- 
' lems are rephrased in terms of difference fields, how the problems are solved there, 

and how the derived results in the given difference field can be reinterpreted as so- 
lutions of the input problem. The algorithms are illustrated with the Mathematica 
\^ ' package Sigma by discovering and proving new harmonic number identities ex- 

C/j , tending those from (Paule and Schneider, 2003). In addition, the newly developed 

C/3 ' package EvaluateMultiSums is introduced that combines the presented tools. 

In this way, large scale summation problems for the evaluation of Feynman diagrams 
in QCD (Quantum ChromoDynamics) can be solved completely automatically. 



> 

^^ ! 1 Introduction 

^, 

l' ' We will elaborate a symbolic summation toolbox based on up-to-date algorithms 

f^*) , in the setting of difference fields. It contains hypergeometric and g'-hypergeometric 

Cn ■ summation, see, e.g., ||Ml[82|52l|58l|52|55][37l|38l|28landi51|^[l6l respec- 

tively, and it can deal with multiple sums covering big parts of (g'-)hypergeometric 
multi-summation ||80l|79l|6i[l5l and (g'-)holonomic sequences ISTI [3U1 l68l l44l . 
. . , This difference field approach started with Karr's theory of J7Z*-fields and his 

r> ' indefinite summation algorithm 1391 l40l : for the continuous analogue of indefinite 



j^ I integration see AMI . In principle, the algorithm solves the telescoping problem in 

a given field generated by indefinite nested sums and products. In this article we 
restrict this input class to nested hypergeometric sum expressions (see Definition^, 
i.e., expressions where the arising products represent hypergeometric sequences, 
and the sums and products occur only as polynomial expressions in the numerators; 
evaluating such expressions produces d' Alembertian sequences lfT2l , ||59l , a subclass 
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of Liouvillian sequences ||36l . We point out that exactly this restriction covers all 
the summation problems that have been relevant in practical problem solving so far 
There we solve the following fundamental problem: given a nested hypergeometric 
sum expression, calculate an alternative expression such that the occurring sums are 
algebraically independent ll76l ; for related work see llTS] |69l W\] [35l |9l . In addition, 
the found representation should be given in terms of sums and products that are 
as simple as possible; for a general framework we refer to [l49l . This is possible 
by representing the sums and products in ITZ* -fields reflecting certain optimality 
properties: We will exploit simplifications taking into account, e.g., the minimal 
nesting depth ||65l|621El|75l|73 or minimal degrees ||69] |72l [131 . 

Besides indefinite summation, we aim at the transformation of a definite multi- 
ple sum to nested hypergeometric sums. As for the special case of hypergeometric 
summation ||58ll82l one looks for a recurrence of such a sum ll62l . If one succeeds, 
one computes all solutions of the found recurrence that are expressible in terms of 
nested hypergeometric sum expressions; for solvers of recurrences in terms of poly- 
nomials and nZ* -fields see |I27l[l2l[36l and |l26l|62l|70l[IIl, respectively. Finally, 
one tries to combine the solutions to an expression that equals the input sum. 

All these algorithms (also for the ^-hypergeometric and mixed case) are available 
in the summation package Sigma |73| and have been used to discover and prove 
demanding identities from combinatorics or related fields, like, e.g., in ll32l 17111561 
|5T1 . A typical example is the sum 

\ « k ^ 



A„(fl) = ^ (1 + a(n - 2k)Siik)) with 5i (k) = ^ - 

jfc=o V'^/ i=\ ' 

which is connected to supercongruences of the Apery numbers. For the treatment 
of the cases a G {1,2,. ..,5} and a > 5 we refer to |!54] and ||45] , respectively. As 
running example we will discover and prove the following identitie^ 

A_i(«)=(n + l)5i(n) + l, (1) 

{n+lf , {n+2{n^ + 3n + 2)Si{n)+3){n + l) 

[n + iy [n + iy 

A_3(«) =(-l)"(55_3(«)(« + 1)^ -65_2.i(n)(n + 1)^) +65i(n)(« + !) + !, (3) 

(4«(« + 2)+3) \2,fi i^ f-{ i^ I 



(4) 



(io(« + i)^i(«)+3)(« + i; 

^ 2« + 3 

where the harmonic sums ll25l l78l are defined by 

A sign(mi)''i 'i-,' signK)" ^ ^ \ rnl /cn 

Snn.......m,{n)= 2^ .|,„^| ■■■ L .\m,\ ' mieZ\{Q}. (5) 

'1=1 'l 't=l 'k 



' For identities (T), J2) we point also to (29); for their indefinite versions see <10t .l llU below. 
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We emphasize that exactly this type of nested hyperg eometri c sums is related to 
summation problems coming from QCD like., e.g., in |48 23 1 or in ifTTl l6l l4l |2T1 . 
More precisely, 2- and 3-loop Feynman integrals with at most one mass and with 
operator insertion can be transformed to multiple sums [g41 depending on a dis- 
crete Mellin parameter n. Then these sums must be simplified in terms of special 
functions ||2l, such as harmonic sums (|5]l, their infinite versions of multiple zeta val- 
ues ||20| and generalizations like 5-sums 1471 |9l and cyclotomic harmonic sums fE\. 
In recent calculations |3 1 also binomial sums as in (|4|i arose. For certain sum classes 
we point to efficient tools like 17811471 . For harder sums such as @| 



/(-3 / k -j+n-3-l+n-q-3-l+n-.,- 

Ill E I E 

j=Ok=QI=Q q=0 



.-.„,.„ ,.„ ,., h (-/+«-,-2)!(-;+«-l)(„-,-r-.-2)(,+.+ l) 

(^_iyJ+i'-i+n-ci-i]^4Si(-j + n-l)-4S,{-j + n-2)-2S,{k)-(Si{-l + n-q-2) (6) 
+ Si{-l + n-q-r-s-3)-2Si{r + s)) + 2Si{s-l)-2Si{r + s)j 

the summation techniques under consideration work successfully and are applied 
automatically in the newly developed package EvaluateMultiSums ll22l . In 
this way, millions of multiple sums l6| |2T| could be treated. In addition, the pre- 
sented packages and algorithms are used by an enhanced version IH O of the 
method of hyperlogarithms ||27| and by new algorithms for the calculation of e- 
expansions 1241 171 utilizing multi-summation and integration methods l|79l [TSl l68l . 
The outUne of this article is as follows. In Section |2] we present the basic mech- 
anism how expressions in terms of indefinite nested sums and products can be 
rephrased in a difference field. For interested readers details are given in Section[3| 
Readers, that are primarily interested in the summation tools and how they can be 
applied with the summation package S igma, can jump directly to Section|4] Finally, 
in Section|5]the new package EvaluateMuliSums is introduced that combines 
all the presented summation methods. It enables one to simplify definite nested sums 
to indefinite nested hypergeometric sums completely automatically. 

2 Indefinite summation: the basic mechanism 



We will work out the basic principles how indefinite summation can be carried out 
in the setting of difference fields. This will be illustrated by the task to simplify 

A_,ia)='£Fik)^t{l-in-2k)S,ik))('l) , (7) 

k=o k=o V''/ 

to be more precise, by the task to solve the following problem. 



Problem T: Telescoping. Given a summand F{k). Find an expression G{k) such 

that 

Gik+l) = G{k)+F{k+l) (8) 

and such that G{k) is not "more complicated" than F{k). 



For our given F{k) in (|7]i we will compute for k>Q the solution 
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G{k) = {{k+l)Si{k) + l){;iy'+c, ceQin). (9) 

Since A_i(a) andG(fl) satisfy both the recurrence A (a + 1) —A{a)+F{a + 1), they 
are equal for a > if they agree at a = 0; this is the case with c = 0. Hence we get 
for A_i (a), and with the same technique for A_2(fl), the simplifications 

A_i(fl) = ((fl + l)5i(fl) + l)Q ', (10) 

, , , {n + lf {a + l){-a + 2n+2{a + l){n + 2)Si{a) + 3) fn\-- ,,,, 
A-2(aj = 7 — Hv^H 7 — ^v^t • ^^^^^ 



(« + 2)2 (« + 2)2 



Note that for the special case a = n this simplifies to ^ and dU, respectively. 

Subsequently, we give more details how this solution G{k) for dH) can be de- 
rived automatically. First observe that the occurring sums can be written in terms of 
indefinite sums and products: for all A: £ N, 

S,{k) = i\ and G)=n^; 
1=1 1=1 

here « is considered as a variable. Now let S^i^ be the shift operator w.rt. k. Then 

using the shift behavior of the summand objects, namely 

J^kn = n, y^k^k+l, ^,Q = f^(«), y,Si{k)^S,{k) + ^^, (12) 

we can write, e.g., F(fc+ 1) again in terms of «, fe, (^,) and5'i(^): 

F{k+\) = y,F{k) = {\-{n-2{k+\)){sm + ^))'j^y{"y- (13) 

We will utilize this property, but instead of working with the summand objects k, 
(^) and Si{k) we will represent the objects by the variables x, b, h, respectively; n 
is also considered as a variable. Here we start with the rational numbers and con- 
struct the rational function fielco F := Q(«)(x)(/7)(/i), i.e., the field of quotients of 
polynomials in the variables n,x,b,li. In this way, ( fTsT l is represented by 

f^{l-{n-2{x+mh + ^)'^b-'e¥. (14) 

^ X + 1 .x + 1 

Finally, we model the shift operator ^j. by a field automorphism C7 : F — > F. 

Definition 1. Let F be a field (resp. ring). A bijective map C7 : F — s> F is called ^eW 
(resp. ring) automorphism if o{a ob) = o{a) o o{b) for all a,/? G F and o g {+, •}. 

Remark. If F is a ring, it follows that C7(0) = 0, (7(1) = 1 and C7(-fl) = -O'(fl) for 
all fl G F. In addition, if F is a field, this implies that a{\/a) = \/a{a) for all a G F*. 

Namely, looking at the shift behavior of the summand objects (fT2] | the automor- 
phism is constructed as follows. We start with the rational function field Q(«) and 
define a : Q(«) -> Q(«) with a{c) = c for all c G Q(«). Next, we extend a to 
Q(n)(x) such that a(x) =x+l. We note that this construction is unique: 



^ Z are the integers, N= {0,1,2,...} ai'e the non-negative integers, and all fields (resp. rings) 
contain the rational numbers Q as a subfield (resp. subring). For a set A we define A* := A \ {0}. 
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Lemma 1. Let F(f) be a rational function field, o :¥ ^¥ be a field automorphism, 
and a,b g¥ with a ^0. Then there is exactly one way how the field automorphism 
is extended to F(f) subject to the relation G{t) = at + b. Namely, for f = Y!i=ofif' G 
ni a{f) = UlMfi){^t + bY- Andforp,q G F^ with q^Q, a(£) = ^. 

As a consequence, by iterative application we extend CJ uniquely from Q(«) to 
Q,{n){x){b){h) subject to the shift relations (compare (fTzt ) 

a{x)=x + \, a{b) = '^b, a{h)^h + -^. (15) 

X + 1 A' + 1 

In summary, we represent the summand F{k+ 1 ) given in ( fT3] l by ( fT4l i in the rational 
function field F := Q{n){x){b){h) together with its field automorphism a : F -^ F 
subject to the shift relations (fTST l. Exactly this construction is called difference field; 
for a general theory see 



Definition 2. A difference field (resp. difference ring) (F, a) is a field (resp. ring) F 
together with a field automorphism (resp. ring automorphism) CJ : F ^ F. Here we 
define the set of constants by consto-F := {c e F|c7(c) = c}. 

Remark. For a difference field (F, CJ) the set const(jF forms a subfield F which is 
also called constant field of (F, cr). Since Q is always kept invariant under a (this is 
a consequence of cj(l) = 1), Q is always contained in const(jF as a subfield. 
We continue with our concrete example. Given the difference field (F, d) in which 
F (^ + 1 ) is represented by (fl4t . we search for a rational function g G F such that 

o{g)^g + f- (16) 

Namely, we activate the algorithm from Section|32]below and calculate the solution 

g = {{x+\)h + \)b-^+c, cgQ(«) (17) 

which rephrased in terms of the summation objects gives the solution (|9]l for (O. 
In a nutshell, the proposed simplification tactic consists of the following steps. 

1 . Construct a difference field in which the summand objects can be rephrased. 

2. Find a solution g of ( fT6b in this difference field (or a suitable extension). 

3. Reformulate g to a solution G{k) of ^ in terms of product-sum expressions. 

The algorithms in the next section deliver tools to attack this problem for the class 
of indefinite nested product-sum expression; for a more formal framework see fT5\ . 

Definition 3. Let K be a field and the variable k be algebraically independent over 
K. An expression is called indefinite nested product-sum expression w.r.t. k iff it can 
be built by k, a finite number of constants from K, the four operations (+, —,•,/), 
and sums and products of the type Lf=//(0 or nf=//(') where / G N and where 
/(/) is an indefinite nested product-sum expression w.rt. / which is free of fc . In 
particular, we require that there is a A G N such that for any integer « > A the 
expression evaluates for k = n (to an element form K) without entering in any pole. 
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In this recursive definition the sums and products can be arbitrarily composed and 
can arise also as polynomial expressions in the denominators. Here we restrict our- 
selves to those expressions that occurred in practical problem solving so far. 

Definition 4. A sequence (/Jj,)h>o G K^ is called hypergeometric if there are a{x) e 
K(x) and / e N such that hu+\/hu = a{u) for all u>l. I.e., /2„ = c n"=/+i ^i^ ~ 1) 
for all M > / for some c G K*. Such a symbolic product (m replaced by a variable k) is 
called hypergeometric product w.rt. k. An expression is called nested hypergeomet- 
ric sum expression w.rt. k if it is an indefinite nested product-sum expression w.rt. k 
(see Def . [3]) such that the arising products are hypergeometric and the arising sums 
and products occur only as polynomial expressions in the numerators. The arising 
sums (with upper bound k) are called nested hypergeometric sums (w.r.t. k). 

E.g., the harmonic sums ^ and their generalizations P7] |8] fall into this class. In 
particular, the right hand sides of ([T)-(|4) are covered. These expressions evaluate 
exactly to the d'Alembertian sequences |[T2l . ||59l . 

Remark 1. In the difference field approach also the ^-hypergeometric and mixed 
case lfT6l can be handled. All what will follow generalizes to this extended setting. 

Subsequently, we will derive a full algorithm that treats the three steps from 
above automatically for the class of nested hypergeometric sum expressions. 



3 Details of the difference field machinery 

We will work out how the Steps 1-3 from above (covering also the more general 
paradigms of creative telescoping and recurrence solving) can be carried out auto- 
matically. As an important consequence we will obtain tools to compactify nested 
hypergeometric sum expressions, i.e., the occurring sums in the derived expression 
are algebraically independent (see also Section in l4.lt . 



3.1 Step 1: From indefinite nested sums and products to TIZ* -fields 

In the previous section the construction of a difference field for a given summand in 
terms of indefinite nested product-sums was as follows. We start with a constant dif- 
ference field (K, a), i.e., a{c) — c for all c G K or equivalently const(jK = K. Then 
we adjoin step by step new variables, say f i , . . . , r^ to K which gives the rational func- 
tion field F := IK(f i ) (^2) ■ ■ -{te) and extend the field automorphism from IK to F sub- 
ject to the shift relations CT(r,) = a,?; or o{ti) = tt + at for some a,- G ¥{t\ ) . . . (r,-i )*. 
Subsequently, we restrict this construction to J7Z*-fields; for a slightly more general 
but rather technical definition of Karr's Hi^-fields see I 



Definition 5. (F, a) as given above is called TIL* -field over IK if constfjF = IK. The 
adjoined elements (f 1 , . . . , f^) are also called generators of the ITZ* -field. 
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E.g., our difference field {Q{n){x){b){h),a) with ([Bl) is a ni:*-field over Q(n). To 
see that the constants are just Q(n), the following result is crucial 



Theorem 1. [Karr's theorem] Let (F,(7) be a difference field, take a rational func- 
tion field¥[t), and extend the automorphism O from F fo F(f ) subject to the relation 
o{t) = at + f for some a G F* and / G F. Then the following holds. 

1) Case a = 1.- const(jF(f ) = constfjF iff there is no g <E¥ with G{g) ~ g + f- 

2) Case f = 0: consto-F(f ) = constg-F iff there is no g e¥*, r>0 with a(g) == a'' g. 

Example 1. Using Theorem[T]we represent the sum (|2) in a m;*-field parsing the 

occurring objects in the following order: —In—lk-i (^) — >5i(A:) — >A_i(a). 

(0) We start with (Q(n), a) setting a(c) = c for all c G Q(«). 

( 1 ) Then we construct the difference field (Q(n) {x),a) subject to the shift relation 
g{x) =x+l. Since there is nog G Q(n) such that a{g) ^ g+l, it follows by Karr's 
Theorem that constaQ(«) W = Q(«). i-^-, (Q(«)(x),(j) is a iTi:*-field over Q{n). 

(2) One can check by an algorithm of Karr lf39l that there is no r > and g G 
Qln){x)* such that a{g) = {'^Y g- Thus for our difference field {Q{n){x){b),a) 
with G{b) = '-^b we have that const(,Q{n){x){b) == consto-Q(n)(x) = Q(«) by The- 
oremlU i.e., (Q{n){x){b),a) is a 77i:*-field over (Q(«). 

(3) Next, we extend the ni:*-field {Q{n){x){b),a) to {Q{n){x){b){h),a) subject to 
the shift relation a{h) ^h + -^. There is nog eQ{n){x){b) with a{g) =8+-^; 
this can be checked by the algorithm given in Section l372l below. Thus the constants 
remain unchanged by Theoremfl] and {Q{n){x){b){h),a) is a J7Z*-field over Q(«). 

(4) Given / in (O, that represents F (A: + 1 ) in ©, we find ([T7]i such that a{g) = 
g + f- In other words, g reflects the shift behavior of A_i(fc) = Lf=i^(0 with 
5^kA-\{k) = A^\{k) +F{k+ 1). Reformulating g in terms of sums and products 
yields (|9]) and choosing c = delivers the identity ( fTOl i. In other words, g (for c~Q) 
can be identified with the sum A_ i {k). This construction will be done more precise 
in Subsection 13. 41 in particular, we refer to Remark 5] 

In general. Theorem [T] yields the following telescoping tactic to represent a given 
indefinite nested product-sum expression (see Definition|3) in terms of a nZ*-field. 
One starts with the constant field (K, a) with o{c)~c for all c G K. Then one parses 
all the summation objects. Suppose one treats in the next step a sum of the form 
Y!i=iF{i) where one can express F{k) in the so far constructed 7TZ*-field (F,(7), 
say F{k+\) can be rephrased by / G F. Then there are two cases: one finds a g G F 
such o{g) = g+ f and one can model the sum ^f^i F{i) with its shift behavior 

k k 

ykt.Fi') = LPiO+F{k+l) (18) 

(=1 1=1 

hy g + c (for some properly chosen c G K). If this fails, one can adjoin a new variable, 
say t, to F and extends the automorphism to a : F(f ) -^ F(f) subject to the shift 
relation cj(f) = t + f. By Theorem[T|the constants remain unchanged, i.e., (F(f ), c) 
is a i7Z*-field over K, and t G F(f ) models accordingly the shift behavior (fTsT l of 
our sum. The product case can be treated similarly; see also Problem RP on pagefTOl 



8 Carsten Schneider 

3.2 Step 2: Solving the telescoping problem in a given TIE* -field 

Karr's algorithm ||39l solves the telescoping problem within a fixed J7Z*-field ex- 
ploiting its recursive nature: it tries to solve the problem for the top most generator 
and reduces the problem to the subfield (i.e., without the top generator). This reduc- 
tion is possible by solving the following more general problem. 



Problem FPLDE: First-order Parameterized Linear Difference Equations. 

Given a 7Ti:*-field (F, a) over K, ao, ai e ¥* and /o, . . . ,/^ G F. 

Find alO cq, . . . , q/ G K and § G F such that tti ct(§) + Oog = cq/o H h cjf^i. 



" The solution setV = {(cq, . . . ,(;v,,g) e K''+' x ¥\aicy {g) + a^g = lt(,t,/;} forms a K-vector 
space of dimension <d + 2 and the algorithm calculates an explicit basis of V . 



Remark. Problem FPLDE contains not only the summation paradigm of telescoping, 
but also of creative (resp. parameterized) telescoping (|40] | for a fixed nZ*-field. 

Subsequently, we sketch a simplified version of Karr's algorithm applied to our 
concrete problem: Given the iTi:*-field (F, a) with F = <^{n){x){b){h) and the shift 
relations ( fTSl l and given the summand (fl4] i. calculate (if possible) g G F such that 
o{g)— g^ f holds. The algorithm is recursive: it treats the top most variable h and 
needs to solve FPLDEs in the smaller 7Ti;*-field (H, a) with H = Q(«)(x)(/7). 

Denominator bounding: Calculate a polynomial q G M.[h]* such that for any g G 
M{h) with (fTSl l we have that gq G H[/2], i.e., q contains the denominators of all the 
solutions as a factor For a general HZ* -field and / such a universal denominator q 
can be calculated; see lf39l l26l |63]| . Then given such a q, it suffices to search for a 
polynomial p G M[h] such that the first order difference equation 

1 a{p)~-p = f (19) 



a{q) q 

holds (which is covered by Problem FPLDE). In our concrete example the algorithm 
outputs that we can choose q= I, i.e., we have to search for a p G W[h] such that 
c{p) —/? = / holds. 

Degree bounding: Calculate b such that for any p G IHI[/2] with ( fT9] l we have that 
deg(p) < b. For a general ni;*-field and / such a b can be calculated; see Il39ll66ll . 
In our concrete example we get b = 2. Hence, any solution p G IHI[/2] of G{p) — p = f 
is of the form p = p2li^ + pih + po and it remains to determine p2,pi,po G H. 
Degree reduction: By coefficient comparison of h^ in 

<7{p2h^ + pih^ + Pq) - iP2h^ + Pih^ + Pq) -- f (20) 

we obtain the constraint <7{p2) — P2 = on p2- Since (H, a) is a nZ*-field, p2 G 
Q(«). Hence we can choose p2 = d where d G Q(«) is (at this point) free to choose. 
Now we move p2li~ = dh^ in ( |20] | to the other side and get the equation 

oipih' +po) - {pxh' +po)^f^d ^'[^ll]^K (21) 
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Note that we accomplished a simplification: the degree of h in the difference equa- 
tion is reduced (with the price to introduce the constant d). Now we repeat this 
degree reduction process. By coefficient comparison of h^ in ( |2TI ) we get the con- 
straint a{pi)— pi = — — llij^Sj] — ^+'^7rr onpi. Again we succeeded in a reduction: 
we have to solve Problem FPLDE in H. Applying the sketched method recursively, 
gives the generic solution d = and pi ~ ^ji +e with e G Q(«). Plugging this solu- 
tion into ( |2TI ) and bringing a{pih) — pih to the right hand side reduce the problem 
to (J{pq) — Pq = 1^^",^) ' +^jTr' ^^^^ ^^^ ^^ decreased the degree of h from 1 to 
0, i.e., we have to solve again Problem FPLDE in H. Recursive application of the 
algorithm calculates the generic solution e = and po = i>^' + c with c G Q(n). 
Putting everything together gives the solution (fTTt . 

The technical details of the sketched algorithm for Problem FPLDE can be found 
in Il65ll26l . More generally, this algorithm can be extended to a method from the 
first-order case to the mth-order case {m e N) as described in ifTOl . Furthermore, 
taking results from ifTTI we obtain a full algorithm that solves the following 



Key problem PLDE: Parameterized Linear Difference Equations. 

Given a ni;*-field (F, a) over K, Oq, . . . , a,„ e F (not all zero) and /o, . . . ,fd £ ^■ 
Find alO cq, . . . , Q £ K and g G F such that 

a,„ d"' (g) + • • • + ao§ = co/o + ■ • ■ + c4d . (22) 



" The solution set V = {{cq,. . .,cj,g) € K''+' xF | <22b holds! forms a K- vector space of 
dimension <m + d+l and the algorithm calculates an explicit basis of V . 



Remark. Problem PLDE covers telescoping (see (fTSIl). creative telescoping (see ( |40| |) 
and recurrence solving (see (l4&t ) for a given nZ*-field. In particular, it is a crucial 
building block for the enhanced summation paradigms given below. Furthermore, it 
allows to deal with holonomic sequences in the setting of difference fields Il68ll73l . 



3.3 Restriction to polynomial TIL* -fields 

We described how the summation objects can be rephrased in a nZ*-field (Step 1) 
and how the telescoping problem, and more generally Problems FPLDE and PLDE 
can be solved there (Step 2). Subsequently, we restrict to polynomial nZ*-fields. 
This will allow us to reformulate the found result completely automatically in terms 
of the given summation objects (Step 3) in Subsection l3.4l 



Definition 6. A i7Z*-field (F,C7) over IK with F = K(x)(/:>i) . . . {pr){si) . . . {se) is 

called polynotnialif <j{x) =x+l, 

• G{pi) = ciipi with a^G K(x)* for all 1 < / < r, and 



I VVlLll LJ, ^ 

a{si) = Si+fi, witlU fi e K{x)[pi,pY\. . . ,p,-,/?,7'][ii, . . . ,s,-i] for all 1 < / < e. 



^ K(jc)[pi,P[ ^ ,. . . ,pr,p,. '] stands for the polynomial Laurent ring in the variables p\,. . .,pr, i.e., 
an element is of the fonn £/,| lAesfUi i,)P'\ ■ ■ -P'r where f^j^ ,■ i € K(-^) and 5 C Z*" is finite. 
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Related to Remark[T]we note that in 11761 a more general definition is used that covers 
also the (7-hypergeometric and mixed case |fT6l . All what will follow generalizes to 
this general setting. Let (F, cj) be a polynomial iTi;*-field over K as in Definition|6] 
and define the ring 

R = K{x)[pi,p^\...,pr,p'^^][su...,Se]. (23) 

Note that for all g e R and k eZwe have that o'^ig) G R- Thus restricting a io R 
gives a ring automorphism (see Definition [T). Therefore {R,(y) is a difference ring 
and the set of constants is the field K. 

Example 2 (SeeEx.\B- {Q{ii){x){b){h),a) is a polynomial ili;* -field overQ(«). In 
particular, we get the difference ring {R, a) with constant field Q(«) for the polyno- 
mial (Laurent) ring 

R = Qin){x)[b,b-'][h]. (24) 

We highlight that polynomial ITZ* -fields cover (up to the alternating sign) all nested 
hypergeometric sums (see Definition lU. This can be seen as follows. 

• Hypergeometric sequences. Consider, e.g., the hypergeometric products 

k k 

Hi{k)^Y\ai{i-l), ... , H„ (k) = ]~[ a,, (/ - 1 ) with aj {x) e K{x) , (25) 

(=/l ;■=/,. 

(Hj{k) ^ for all A: > 0) where K = Q or K = Q(ni, . . . ,n„) is a rational function 
field. Then there is an algorithm ll69l Theorem 6.10] based on Theorem[T|that solves 



Problem RP: Represent Products. 

Given a iTi;*-field (K(x),(t) over K with a(x) =x+ 1; ai,... ,«,, e K(x)*. 
Find a JTi:*-field (F, o) over K with F = K{x){pi) . . . (pr) and a{pi)/pi e K{x) 
for 1 < ! < r together with gj G 'K{x)\p\,p^^ ,. . . ,pr,p^^]* andbj e {-1,1} for 
1 < 7 < V such that o{gj) ~ bjUjgj. 



Namely, given ai,...,a,, S K{x), let (F,(7) with R := K{x)[pi,p]^\...,pr,p^^] 
together with gj G R and bj G { — 1,1} for 1 < 7 < v be the output of Problem RP. 
If bj = 1 for all 1 < 7 < V, the products Hj{k) in ( |25] | can be expressed with Cjgj for 
appropriate Cj G K* in the polynomial 7TZ* -field (F, a). 

Otherwise, construct the difference ring (F[m],c7) with a{m) = — m, m' = 1 and 
constaF[m] = K; see ^. Here m models (-1)'^ with S^k{~l)'' = -(-I)'''- Then 
the Hj{k) in (l25T l are rephrased with Cjm^^^^j'^'^gj for appropriate Cj G K*. 

• Indefinite nested sums. Take an expression in terms of nested hypergeometric 
sums, i.e., the sums do not occur in a denominator Moreover, suppose that all the 
arising hypergeometric products can be expressed in a polynomial 7TZ*-field. Thus 
it remains to deal only with summation signs and to extend the given polynomial 
J7Z*-field using Theorem [T]l. Suppose that during this construction it was so far 
possible to obtain a polynomial nZ*-field, say it is of the form (F, a) as given in 
Definition |6] with the difference ring [R, a) with R as in ( |23] |. and let f ^ R be the 
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summand of the next sum under consideration. Then there are two cases. If we fail 
to find ag e F such that a{g) =g + f then we can construct the i7i;*-field (F(f ), a) 
with a{t) = f +/ by Theorem[T] In particular, this i7Z*-field is again polynomial. 
Otherwise, if we find such a § gF with a{g)—g = fGR, we can apply the following 
result; the proof is a slight extension of the one given in ll76l Thm. 2.7]. 

Theorem 2. Let (F, <j) be a polynomial TIL* -field over IK and consider the differ- 
ence ring {R, (j) as above. Let g €¥. If o{g) ^ g £ R then g (z R. 

Example 3 (Cont. Example^. For / G /? with (fT4l i it follows that any solution ^ G F 
with G{g) = g -\- f is in R. Indeed, we calculated ( fTTb . 

Thus we always have g £ R. Asa consequence we can express the sum over / with 
g + c for some properly chosen c G K (see Remark|4]l. Hence by iterative application 
of the above construction we never enter in the case that sums occur in the denomi- 
nators. Consequently, a nested hypergeometric sum expression can be rephrased in 
a polynomial HZ* -field up to the following technical aspect. 



Remark 2. If the hypergeometric products (|25> cannot be expressed in a nZ*-field 
solving Problem RP, one needs in addition the alternating sign (—1)*^ in the setting 
of difference rings; note that here we cannot work anymore with fields, since zero 
divisors pop up: ( 1 + ( — 1 )*') ( 1 — ( — 1 )*^) = 0. For simplicity, we restrict ourselves to 
polynomial nZ*-fields; the described techniques and algorithms in this article can 
be extended for the more technical case allowing also (— 1)*^; see 



3.4 Step 3: Evaluating elements from a TIE* -field to sequences 

Let (F, (j) be a polynomial FIE* -Held over K as in Definition|6]and define R by ( 
In this section we make the step precise how elements from R can be reformulated 
as a nested hypergeometric sum expression. I.e., how such an element f £ R can 
be mapped via an evaluation function ev : i? x N — ?► K to a sequence (ev(/, /)),>o by 
using an explicitly given nested hypergeometric sum expression w.rt. a variable k. 

Before we start with a concrete example, we emphasize that this map ev should 
respect the ring structure R and the ring automorphism a as follows. 

Definition 7. A map ev ; i? x N ^- K for a difference ring {R,(y) with constant field 
K is called evaluation function if the following evaluation properties hold: For all 
c G K and all / > we have that ev(c, /) = c, for all f,g £ R there is a 5 > with 

V; > 5 : ev(/^,/) = ev(/,/)ev(g,0, (26) 

V/ > 5 : ev{f + g,i) = ev(/,0 +ev(g,0; (27) 

and for all / G /? and y G Z there is a 5 > with 

V/>5:ev((J^(/),/)=ev(/,/ + ;). (28) 
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Example 4 (See Ex. [7]). Take the polynomial ni:*-field (Q(«)(x)(fe)(/i),a) with 
(fTSI l and consider {R, a) with (l24l i and constant field Q(n). We construct an evalua- 
tion map ev : 7? X N — > Q(n) as follows. For / G ^ and A: G N, ew{f,k) is calculated 
by evaluating an explicitly given nested hypergeometric sum expression. 
The construction is performed iteratively following the tower of extensions in R. 

(1) We define ev : Q(«)(x) x N ^ Q(«) as follows. For f G Q(n)(x) with p,q G 
Q{n)[x] andgcd(/5,^) = 1, 

e,(^*)Jil '"W^" (29, 

^ [0 ifq{k)=0 (pole case); 

here p{k),q{k) with A: G N denotes the evaluation of the polynomials atx = k. Note 
that the properties in Definition |7] are satisfied for 5 chosen sufficiently large: that 
is the case when one does not enter in the pole case in ( |29] | for elements f,g^Ras 
given in Definition|2l see also Example|5] 

(2) Next, we extend ev from Q{n){x) to Q{n){x)[b,b'^]. We set 

ev{b,k)=cifl "^^,^' , /i gN,ci gQ(«)* (30) 

and prolong the ring structure as follows: for / = Yfj=„fjbJ G Q{n){x)[b,b^^] with 
/;■ G Q(n)(x) and u,d G Z we define ev{f,k) — Yf''^^ev{fj,k)ev{b,k)J; this implies 
that ( l26b and (|27] | hold for some 5 G N sufficiently large; see Example |5] Note that 
for any choice of /i G N and c G Q(n)* also ( l28T l is valid. Since we want to model 
Q = nf^j n+t^, a natural choice is h = 1, ci = 1. 

(3) Finally, we set 

'' 1 
ey{h,k)= y-+du AiGN,^iGQ(«). (31) 

■J ' 

Again properties (l26t and (l27T i hold for some 5 G N (see Example |5]l if we extend 
ev as follows: for / = Tf-^ofjhJ G R with // G Q{n){x)[b,b^^ we set ev(/,A;) = 
Yfi=Q^'^{fijk)ew{h^ky . In addition, property (l2Ft holds for any choice of Ai G N 
and d\ G Q(«). Since we want to model ^i {k), we take, e.g., Ai = 1 and d\ = 0. 
For instance, for / G /? as in ( fT4l) the evaluation is given by the nested hypergeomet- 
ric sum expression (cf. ( fTJt ) 

ev(/,^)=(l-(n-2(^+l))(5,(fc) + ^))^Q ; (32) 

the usage of Si {k), ('^) is just pretty printing and stands for (l30t , d3Ti . respectively. 

Besides the function ev we aim at the calculation of the bounds 8 in Definition|2l 

Example 5. For the evaluation function ev : 7J x N — )► (Q(«) from Example |4] the 
bounds can be extracted by the computable map P : R ^-N defined as follows. For 
/ G ^ let (i G N be minimal such that for all k > d the calculation of e\{f,k) 
does not enter in the pole case in (|29] l. More precisely, write / in the form 
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/ = Y.iez.jenfi.jb'^'^'' with finitely many fjj G Q{n){x) being non-zero, and choose 
d €N such that for all /, j and all ^ > rf the denominator of fij has no pole 3.tx = k. 
This d can be calculated explicitly for any f E R and defines the function j3 with 
J3(/) := d. Now let f,g e R. Then for 5 := max()3(/),j3(g)) we have that 
and JTU . In addition, for all 7 G Z choose 5 := j3 (/) + max(0, —j) and we get ( 

This example motivates the following definition ll76l . 

Definition 8. Let {R,(y) be a difference ring with constant field K and consider an 
evaluation function ev : /? x N ^ K. j3 : /? — > N is called bounding function of ev if 
for all f,geRwe can take 8 := max()3 (/), jS (g)) such that (HH and (|27]) holds, and 
for all / G /? and ; G Z we can take 5 := j3 (/) + max(0, -;) such that (|28ll holds. 

The concrete construction above carries over to the general case. Let (F, a) be a 
polynomial J7Z* -field over IK as given in Definition |6] and define R by (|23] |. Then 
one obtains an evaluation map ev : /? x N — !> IK in terms of explicitly given nested 
hypergeometric sum expressions as follows. We start with ev : K(x) x N — > K de- 
fined by ( [29] l for - G K(x) with p,q G K[x] and gcd{p,q) — 1. Next we define how 
the map acts on the pf. 

k 

ev(p;,^) = c;P]«K;-i), (1 <'•<'■); (33) 



here we are free to choose c, G IK*, and /,■ is chosen such that the numerator and 
denominator of a, evaluated at j is non-zero for all j > /,. Then this map is ex- 
tended to ev : ^0 x N — > K with ^0 := IK(x)[pi,pj"',... ,79,.,p,7^] as follows. For 

/ == I(i,,...,;,.)GZ' /(/i i,.)P'i ■ ■ -P'r e Ro with /(,-, ;^) G K{x) we set 



ev(/,A:)= Y. ev(/(,•,,...,■.,,fc)ev(/7l,^)'l...ev(p„^)'^ 

(n ir)€Z'- 

Finally, we extend iteratively this map from ^0 to R := Re- Suppose that we are 
given already the map for R, = R[si,... ,i,-i] with 1 < / < e. Then we define 

k 

ev(.,-,^)- ^ev(/,-,;-l)+flf,- (34) 

i=A, 

where dt G IK can be arbitrarily chosen, and A, G N is sufficiently large in the fol- 
lowing sense: it is larger than the lower bounds of the arising sums and products 
of the explicitly given nested hypergeometric sum expression for e\{fi,j — 1) and 
such that during the evaluation one never enters in the pole case in ( |29] l. In a nut- 
shell, the underlying expression for ev{si,k) (for symbolic k) can be written as an 
indefinite nested sum without entering poles that are captured via (|29] l. Finally, we 
extend this construction to ^,_i [si]: for / = L/=o//'^/ ^ ^i-i i^'] ^i'^^ fj ^ ^'-i ^^ 
define sv{f,k) ~ L}=o^^(//'^)^^('*i'^)''- To this end, by iteration on / (1 < / < e) 
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we obtain ev : RxN ^K with R := R^ that satisfies the properties in Definition |7] 
and which is exphcitly given in terms of nested hypergeometric sum expressions. 

Remark. Note that the products in ( l33T l and sums in ( |34] | are just nested hypergeo- 
metric sum expressions w.r.t. k; see DefinitionlH 

Moreover, we can define explicitly a bounding function j8 : /? — s> N that produces the 
required bounds 5 in Definition|2]following the construction in Example|5] In short, 
consider f E R as a polynomial in the variables pi,Si and take all its coefficients 
from K(x). Then )3(/) is the minimal value d gN such that the evaluation of the 
coefficients does not enter in the pole case of (l29] i: note that the positive integer 
roots of the denominators can be detected if IK is computable (in particular, if one 
can factorize polynomials over K). Summarizing, we get Lemma|2l cf. 1761 . We can 
summarize this construction as follows. 

Lemma 2. Let (F,(7) be a polynomial Till* -field over K with R as in ( I23l l. Then 
ev : 7? X N — > K defined above in terms of nested hypergeometric sum expressions 
is an evaluation fimction, and j3 : 7? — > N given above is a corresponding bounding 
fi^nction ofev. IfM. is computable, such functions can be calculated explicitly. 

Remark 3. Solving Problem PLDE. Take a polynomial nZ*-field (F, a) as in Def- 
inition |6] and define R by ( |23] ). Let ev : /? x N ^- K be an evaluation function with 
a bounding function j3 : /? — > N as constructed above by means of hypergeometric 
sum expressions. Let a,,/, e R and take g GR such that equation (l22l i holds. Now 
calculate 5 := max(j3(o!o), ... ,j3(a,„),j3(/o), ... ,j3(/j),j3(g)). Then for any fc > 5, 

ev(a„,,^)ev(§,^ + m)H Vc\{aQ,k)&\{g,k) ^ CQev{fQ,k) + ■ ■ ■ + c,tev{fd ,k). 

(35) 
Hence a solution g E R of ( l22l i produces a solution of (|35] | in terms of nested hy- 
pergeometric sum expressions. In particular, a lower bound 5 for its validity can be 
computed (if K is computable). We emphasize that this property is crucial for the 
automatic execution of the summation paradigms given in Section|4]below. 

Remark 4. Simultaneous construction of a TIE* -field and its evaluation function. 
In order to model the summation problem accordingly (see, e.g.. Problem EAR on 
pagefTTIi the construction of the TIE* -Held (Step 1) and the evaluation function with 
its bounding function should be performed simultaneously. Here the choice of the 
lower bounds and constants in ( |33] | and ( |34] | are adjusted such that the evaluation of 
the introduced products and sums agrees with the objects of the input expression; 
for a typical execution see Example|4] In particular the evaluation function is crucial 
if a sum can be represented in the already given ni;*-field by telescoping. This 
is, e.g., the case in (5) of Example [T] We succeeded in representing the summand 
F{k) in ( fT3l ) by the element f G Rasin (fT4l i. Namely, using the evaluation function 
from Example m we have ( [32] i for all A: > 0. Then we calculate the solution ( [TT] ) of 
a{g) =g-\-f. In particular, we obtain ev(g, k) = {{k+l )Si {k)-\-l) ('^) " + c with 
c e Q(n) such that for all A: > we have that ev{g,k+l) = ev(g,A:) +ev(/,A:+ 1); 
see Remark[3] Now we follow the same arguments as in the beginning of Section|2] 
A_i(a) and ev{g,a) satisfy the same recurrence A (a + 1) = A(fl) +ev(/,fl + 1) for 
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all fl > and thus ev{g,a) = A_i (a) when choosing c = 0. In this way, we represent 
A_i (fl) precisely with g (c ^0) in the given i7£*-field and its evaluation function. 



3.5 Crucial property: algebraic independence of sequences 

Take the elements from a polynomial i7Z*-field, rephrase them as nested hyperge- 
ometric sum expressions, and evaluate the derived objects to sequences. The main 
result of this subsection is that the sequences of the generators of the nZ*-field 
are algebraically independent over each other To see this, let (F, a) be a polyno- 
mial ni:*-field over K with F = K{x){pi). . . ipr){si). . . (s,) and define R by ^. 
Moreover, take an evaluation function ev : 7? x N — > IK following the construction of 
the previous subsection. Then we can define the map t : R ^ IK^ with 

T(/) = (ev(/,^)),>o - (ev(/,0),ev(/,l),ev(/,2),...). (36) 

Now we explore the connection between the difference ring {R,(y) and the set of 
sequences t{R) = {t(/)|/ G R}. First, we introduce the following notions. 

Definition 9. Let {Ri , Gi ) and {R2,(y2) be difference rings. 

• If /?! is a subring of R2 and ai(/) = <72(/) for all f £ Ri then {Ri,a) is called 
sub-difference ring of {R2,<j). 

• AmapT:/?i ^>/?2 is called ring homomorphism if T(/g) = t(/)t(^) and t(/ + 
g) = t(/) + T{g) for all f,g G Ri- If T is in addition injective (resp. bijective), T 
is called ring embedding (resp. ring isomorphism). Note: if T is an isomorphism, 
the rings Ri and R2 are the same up to renaming of the elements with T. 

• A map T : Ri -^ R2 is called difference ring homomorphism (resp. embedding/ 
isomorphism) if it is a ring homomorphism (resp. embedding/isomorphism) and 
for all « £ Z, / g7?i we have t(c7"(/)) = Cj (t(/)). Note; if Tis an isomorphism, 
(/?i , a) and {R2 , cr) are the same up to renaming of the elements by T. 

With component-wise addition and multiplication of the elements from K^ we 
obtain a commutative ring where the multiplicative unit is 1 = (1,1,1,...); the field 
K can be naturally embedded by mapping A: G IK to k = (A:,fc,A:, . . . ). 

Example 6. Let K(x) be a rational function field, take the evaluation function 
ev : K{x) x N ^ K(x) defined by ^, and define t : IK(jc) -> K^ as in (O. Now 
define the set 

F := T(K(x)) = {(ev(/,fc)),>o|/ G K(x)}. 

Observe that f is a subring of K^. However, it is not a field. E.g., if we multiply 
ev(x,A:)}|(.>o with ev(l/x,A:)}|(.>o, we obtain (0,1,1,1,...) which is not the unit 1. 
But, we can turn it to a field by identifying two sequences if they agree from a 
certain point on. Then the inverse of ev(x,A:)}^>o is ev(l/jc,A;))/.>o. More generally, 
for/ G K(x)* we get ev(/,^))^>oev(l//,^))^>o = 1. 
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To be more precise, we follow the construction from 11581 Sec. 8.2]: We define an 
equivalence relation ^ on K^ by (fl„)„>o ~ (^«)n>o if there exists a 5 > such 
that a„ = bn for all n > 5. The equivalence classes form a ring which is denoted 
by ^(K); the elements of ^(K) (also called germs) will be denoted, as above, by 
sequence notation. Finally, define the shift operator y : S(K) -^ S{K) with 

^((fl0,fll,fl2,...)) = (ai,fl2,fl3,---)- 

In this ring the shift is invertible with o5^^'((fli,a2, •••)) = (0,fli,fl2,fl3, . .•) = 
(flo, fli , fl2, . . . ) . It is immediate that .5^ is a ring automorphism and thus (5(K), ^) 
is a difference ring. In short, we call this difference ring also ring of sequences. 

Example 7 (Cont. Example^. Consider our subring F of ^(K). Restricting S^ Xo F 
gives a bijective map and thus it is again a ring automorphism. Even more, since F 
is a field, it is a field automorphism, and (F, 5^) is a difference field. In particular, 
(F,^) is a sub-difference ring of [S{K),y). 

More generally, consider the map T : /? — > ^(IK) as in (|36] |. Since ev : 7? x N — ^ IK 
has the properties as in Definition|7] it follows that for all f,g(ERwe have T{fg) = 
t(/) T{g) and T{f + g) = t(/) + T{g). Hence T is a ring homomorphism. Moreover, 
for aWfeR and all n e Z, 

^"((ev(/,fc)),>o» = (ev(/,fc + «)),>o) = (ev((j"(/),fc)),>o). 

Thus T is a difference ring homomorphism between {R,(y) and (5'(K),o5^). Since 
t{R) is a subring of S{K) and J^ restricted to t{R) is a ring automorphism, 
{t{R),S') is a difference ring, and it is a sub-difference ring of {S(K),S'). 

Examples. Take the polynomial ITZ* -field {Q{n){x){b){h),a) overQ(«) with O 
and ( |24| ). and let ev : /? x N — > Q(«) be the evaluation function from Example |4l 
define T : R ^ S{Q{n)) with (|36] |. Then T is a difference ring homomorphism. In 
particular, {t{R),S') is a difference ring and a sub-difference ring of (5(Q(n)),o5^). 

Now we can state the crucial property proven in ll76l : our map ([36] l is injective. 

Theorem 3. Let (F, cr) be a polynomial TIL* -field over K., define R bv i23i , and take 
an evaluation fimction ev : /? x N — )■ K as given in Lemma |2] Then T. R ^f 5(K) 
with ( 1361 ) is a difference ring embedding. 

Example 9. (<Q(n)(x),(7) and {T{Q{n){x)),a) are isomorphic. In addition, the rings 
{R,a) and (t(7?),(7) with R := Q{n){x)[b,b-^][h] are isomorphic. Thus t{R) = 
T{Q{n){x))[T{b),T{b^^)][x{h)] is a polynomial ring and there are no algebraic rela- 
tions among the sequences T{b), T(fe^'), t(/i) with coefficients from t(Q(«)(x)). 

In general, the difference rings {R,Ct) and {t{R),J^) are isomorphic: they are the 
same up to renaming of the elements by T. In particular, we get the polynomial ring 

T{R) = T{K{x)Mpi),T{Pl'),---,^iPr)MP7')M^l),---,<^e)] (37) 

with coefficients form the field t(]K(x)). I.e., there are no algebraic relations among 
the sequences t(/7,), '^{Pi) and T(i,) with coefficients from t(K(x)). 
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4 The symbolic summation toolbox of Sigma 

In the following we will give an overview of the symbolic summation toolbox that 
is available in the Mathematica package Sigma 1731 . Here we focus on nested hy- 
pergeometric sum expressions (w.r.t. k) as given in Definition |4l the products are 
hypergeometric expressions (for more general classes see Remark[TJ and the sums 
and products do not arise in the denominators. 

Concerning indefinite summation it is shown how a nested hypergeometric sum ex- 
pression can be compactified such that the arising sums are algebraically indepen- 
dent and such that the sums are simplified concerning certain optimality criteria. 

Concerning definite summation the package Sigma provides the following toolkit. 
In Section l4~2l it is worked out how a recurrence can be computed with creative tele- 
scoping for a definite sum over a nested hypergeometric sum expression. Moreover, 
in Section l43] it is elaborated how such a recurrence can be solved in terms of nested 
hypergeometric sum expressions which evaluate to d'Alembertian sequences. Usu- 
ally the derived solutions are highly nested, and thus indefinite summation is heavily 
needed. Finally, given sufficiently many solutions their combination gives an alter- 
native representation of the definite input sum. Summarizing, the following "sum- 
mation spiral" is applied | 




Remark. We give details how these summation paradigms are solved in the setting 
of polynomial ITZ* -fields introduced in Section [3] These technical parts marked 
with * can be ignored if one is mostly interested in applying the summation tools. 



4.1 Simplification of nested hypergeometric sum expressions 

All of the simplification strategies of Sigma solve the following basic problem. 



Problem EAR: Elimination of algebraic relations. 

Given a nested hypergeometric sum expression F{k). 

Find a nested hypergeometric sum expression F{k) and A G N such that F{k) = 

F{k) for all k>X and such that the occurring sums are algebraically independent. 
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The following solution relies on Section|3]utilizing ideas from 1139117011761 . 
Solution*. Compute a polynomial nr*-fielcO (F, c) over IK as in Definitionl23lwith 
R defined as in ^ together with an evaluation function ev : 7? x N — > IK in which one 
obtains an explicit / G /? with A G N such that ev(/, k) = F{k) for all ^ > A. Output 
the nested hypergeometric sum expression F{k) that encodes the evaluation ev(/, k). 
Concerning the algebraic independence note that the sub-difference ring (l37t of 
the ring of sequences {S{K),S') forms a polynomial ring; here the difference ring 
embedding T is defined by ( l36l l. The sequences given by the objects occurring in 
F{k) are just the the generators of the polynomial ring d37t . 

Remark. K is the smallest field that contains the values of F{r) for all r G N with 
r > X. Here extra parameters are treated as variables. However, in most examples 
these parameters are assumed to be integer valued within a certain range. In such 
cases it might be necessary to adjust the summation bounds accordingly. 

A typical instance of Problem EAR is the simplification of the sum (|7]i: Using the 
above technologies, see Examples[T]and|4]for further details, we can reduce the sum 
A_i (fl) in terms of (") and ^i (a) and obtain the simplification (ITOl i. After loading in 

in[i]:- << Sigma. m 

Sigma — A summation package by Carsten Schneider© RISC 

this task can be accomplished with the function call 

in[2]:=SigmaReduce[^{l-(n-2k)S|l,k])(J) ,a] 

"=" /'n\ "' 
Outpi- ((a+l)S|l,a| + l)( j 

Note that S[mi,... , mic,n] stands for the harmonic sums (|5}. More generally, one gets 
reduced representations for nested hypergeometric sum expressions such as 

k=l k=l k=l k=l k=l 



I Cl^Y-L (V,') = ii + T^a(2a+l)=(4a + 45)(-) = 



Out[3]= 

Simplification witli improved difference field theory 

The solution of Problem EAR is obtained by calculating a set of algebraic inde- 
pendent sums (the generators of the 7TZ*-field) in which the occurring sums of the 
input expression can be rephrased. In order to guarantee that the output expression 
consists of sums and products that are simpler (or at least not more complicated) 
than the input expression, the generators of the iTi^* -field must be constructed such 
that certain optimality criteria are fulfilled. In short, we refine Problem EAR using 
improved iTZ*-difference field theory and enhanced algorithms for Problem T. The 
most useful features of SigmaReduce can be summarized as follows. 
• Atomic representation. By default all sums are split into atomic parts (using partial 
fraction decomposition) and an algebraic independent representation of the arising 
sums and products is calculated. In addition, Sigma outputs sums such that the 
denominators have minimal degrees w.rt. the summation index (i.e., if possible, the 
denominator w.rt. the summation index is linear). A typical example is 



"^ As observed in Remark |2]one miglit need in addition the alternating sign to represent all hyper- 
geometric products. The underlying solution works analogously by adapted algorithms. 
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o- T. a ,T- f -2 + k (l-4k-2k^)S[l,k] (l-4k-2k2)S[3,k] n, , 

,n[4,:= S.gmaReduce[i: ( io(TTF) + 10(l + k^)(2 + 2k + k^) + 5(l + k^)(2 + 2k + k^) )'"' 

10(a2 + 2a + 2) ^ ' ^ 5(a2+2a + 2) ^ ' J = ^ ' %%j k^ 

This feature relies on algorithms refining those given in ll72l : for the special case 
of rational sums see, e.g, lfT4ll52l . By default this refinement is activated; it can be 
switched off by using the option SimpleSumRepresentation->False. 

• By default the following fundamental problem is solved: 



Problem DOS: Depth Optimal Summation. Given a nested hypergeometric 
sum expression. Find an alternative representation of a nested hypergeometric 
sum expression whose nesting depth is minimal. Moreover, each sum cannot be 
expressed by a nested hypergeometric sum expression with lower depth. 



The solution to this problem is possible by the enhanced difference field theory of 
depth-optimal nZ*-fields and the underlying telescoping algorithms; see Il74ll75]| . 
E.g., we can flatten the harmonic sum 53. 2.1(a) of depth 3 to sums of depth< 2: 



in[51:= SigmaReduce[^ ^ L ^ L U 

i=I '" J=I J k=l " 



Out[5]= 



I ^ i f +( t ^)( I ^ ii f )- 1 ^( I ^)( I f 

i-|^ = l 1 i2 1 -'-I ■'- -'-1 ■'- -'-2 1 ^2 ^2 



This depth-optimal i7Z*-field theory yields various structural theorems ITTTI . i.e., 
gives a priori certain properties how the telescoping solution looks like. In partic- 
ular, this leads to very efficient algorithms (for telescoping but also for creative 
telescoping and recurrence solving given below) where we could work with more 
than 500 sums in a depth-optimal nZ*-field. The naive (and usually less efficient) 
lTi;*-field approach is used with the option SimplifyByExt-H-None. 

Example 10. For the 2186 harmonic sums (|5} with weight ^JLj |m,| < 7 all algebraic 
relations are determined filpj . More precisely, using their quasi-shuffle algebra the 
sums could be reduced by the HarmonicSums package IT] to 507 basis sums. 
Then using the algorithms above we showed that they are algebraic independent. 

• Reducing the number of objects and the degrees in the summand. The depth- 
optimal representation can be refined further as follows. 

Given a nested hypergeometric sum expression, find an alternative sum represen- 
tation such that for the outermost summands the number of occurring objects is 
as small as possible (more precisely, concerning a given tower of a iTZ*-field the 
smallest subfield is searched in which the summand can be represented); see [|651 . 

E.g., in the following example we can eliminate ^i {k) from the summand: 

in[6]:= SigmaReduce[ ^ {-l)''Si (k)^ ( " ) ,a, SimplifyByExt -y DepthNumber] 

Ou,[6H -(a-n)(n=Si(a)^+2nSi(a) + 2)<::i^-^_i f^ '^-^^^('^] 

Furthermore, one can calculate representations such that the degrees (w.rt. the top 
extension of a J7Z* -field) in the numerators and denominators of the summands are 
minimal 1721 . For algorithms dealing with the product case we point to 16911131 . 



20 Carsten Schneider 

4.2 Finding recurrence relations for definite sums 

Given a sum, sajo A(n) = L<;=o '""'" F„{k) where F„{k) is a nested hyperge- 

ometric sum expression depending on a discrete parameter n, find polynomials 
co(n), . . . ,Cii{n) (not all zero) and an expression li{n) in terms of sums that are sim- 
pler (see below) than the sum A(«) such that the following hnear recurrence holds: 

co{n)A{n) -\ \-c,i{n)A{n + d) = h{n). (38) 

We treat this problem by the following variation of creative telescoping ll82ll . 



Problem CT: Creative telescoping (general paradigm). 

Given d gN and F„{k) such that F„+i{k) with / G N (0 < / < t/) can be written 
as nested hypergeometric sum expression. Fitid X G N, cq, . . . ,q e K(«), not all 
zero, and G{k) such that for all A: > A we have 

cqF„ {k) +ciF„+i{k) + --- + cjF„+d{k) = 5^kG{k) - G{k) (39) 

and such that the summands of the occurring sums in G(k) are simpler (depending 
on the chosen strategy, see below) than Fn+i{k)\ if this is not possible, return L. 



The following solution relies on 1391 ITOl l76l . 
Solution*. We consider the parameter « as variable. Compute an "appropriate" poly- 
nomial iTi;*-field (F, a) over the constant field K(n) as in Definitionl23]with R de- 
fined as in ^ together with an evaluation function ev : 7? x N — ^ K in which one ob- 
tains explicitly /o ,..., /^ €7? with A' eNsuchthatev(/,-,A:) ^F„+i{k) forallA:> A'; 
again we point to Footnote|4] Compute, if possible, a solution cq, . . . ,c^ G ]K(n) (not 
all zero) and g GR {ox an extension of [R, a) with an extended evaluation function 
ev) such that 

CQh + --- + cdfd^a{g)-g (40) 

holds; in addition we require that in g the summands of the occurring sum exten- 
sions are simpler (depending on the chosen strategy, see below) than each of the 
given/,. If there is not such a solution, return _L. Otherwise extract a nested hyper- 
geometric sum expression G{k) such that G{k) = t\{g,k) and compute A G N such 
that ( l39l ) holds for all A > A; see Remark [5] with m = 1, ai = l,ao = -1. Then 
return (cq, . . . ,cj,G{k)) and A. 

Application. Usually, one loops over li = 0, 1 , . . . until a solution for (|39] l is found; 
for termination issues see Remark 1 on page |26] Then summing (l39b over a valid 
range, e.g., from A to a, gives 

co(«)f F„(A) + ...+Q(«)f F„+^(A)=G(fl + l)-G(A) (41) 

k=X k=X 

where by construction the summands of the arising sums in h{a) :~ G{a + 1) — 
G(A) are simpler than Fn+i{k). This implies that also the arising sums h{a) are 
simpler than LiL;t F,+i{k). Note that so far n is considered as an indeterminate. We 



' L(mi ,... ,m„,n),u>0, stands for a linear combination of the m, and n with integer coefficients. 
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remark that in many applications n itself is an integer valued parameter and extra 
caution is necessary to avoid poles when summing up (|39] l. Finally, when setting 
a = L{mi , . . . ,m„,«) in W\\ (if a = oo, a limit has to be performed) and taking care 
of missing summands yields dJSl ) for A(«); see ExamplefTTIfor details. 

Proof certificate. The correctness of (l3Ft for a given sum A(«) is usually hard to 
prove. However, given the proof certificate (cq,... ,Cij,G{k)) it can be easily verified 
that (|39] | holds within the required summation range. Then summing this equation 
over this range yields the verified result i 



With Sigma one can calculate forA_3(«) = SUM [n] a recurrencelj as follows: 

in[7];= mySum = ^^(1 -3(n-2k)Si(k)) ( 



in[8]:- rec — GenerateRecurrencefmySum.il] [[1]] 

Out[8]= (n + 2)''(n + 3)^SUM[n] + (n+l)^(n + 3)^(2n + 5)SUM[n + l] + (n+l)^(n + 2)^SUM[n + 2] 
== (20n^ + 138n^ + 311n + 229)(n+l)^ + 6(n + 2)^(n + 3)(2n + 5)(n+l)^Si(n) 

The essential calculation steps are given in the following example. 



,-3 



Example 11. * Take A_^{n) = YJl^QF„{k) with F„{k) = (1 - 3(« - 2A:)5i(yt))(^')" 
We calculate a recurrence for A_3(n) in n by the techniques described above. First, 
we search for a solution of ( |39] l with d = Q (which amounts to telescoping). I.e., 
we construct the polynomial iTZ*-field (Q(n)(x)(Z7)(/j), a) and evaluation function 
ev:7?xN^Q(n) with R := 'Q{n){x)[b,b-^][h] as in Example g] There we take 
/o = (1 - 3(2« -x)h)b-^ e R such that ev(/o) = F„{k) for all ;t > 0. Unfortunately, 
our telescoping algorithm fails to find a g E R such that (y{g) — g = fo holds. So we 
try to find a solution of (|39] | with d = 1. Since ("]^ ) = ^"^' j (^'), we can rephrase 

F„+i{k) by /i ^ {l-3{2n-x)h)^^i^^^b-\ i.e., ev(/i,fc) = F„+i(^) for all k>0. 
Then we activate the algorithm for Problem FPLDE and search for cq , c i G Q (« ) (not 
both zero) and g G R such that (|40] | holds with d = I. Again there is no solution. 

We continue our search and take /2 = (1 — 3(2n —x)h) ^"^^^'.S"^l'^^' b^^ with 

ev{f2,k) = F„+2{k) and look for cq,ci,C2 G Q(n) (not all zero) and g G R such 
that ( |40] | holds with d = 2. This time our algorithm for Problem FPLDE outputs 
Co = (n + 2)4(n + 3)2, ci = (« + 1)3(« + 3)2(2n + 5), C2 = in+l)^n+2f and 
8 = {p\{n,x)+ p2{n,x)h)b^^ for polynomials /5i(n,.x),/72(n,-^) &Q[n,x]. Hence we 
get ^ with G{k) = ev{g,k) = {p{n,k) + p2{n,k)Si{kj){"^^^ . We emphasize that 
the correctness of (|39] l for the given solution for all k with Q <k <n can be verified 
easily. Finally, summing ( [39] l over ^ from to « one gets 

co(n)A3(n)+ci(«)(A3(n + l)-F„+i(« + l)) 

+ C2(«)(A3(« + 2)-F„+2(n + l)-F„+2(« + 2)) = G(« + l)-G(«); 

moving the F„+,(n + y) terms to the right hand side gives the recurrence Out(8]. 



* For a rigorous verification the proof certificate (t'o, . . . , q, G(^)) of | |4U with d = 2 is returned 
with the function call CreativeTelescoping [mySum, n] . 
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Note that creative telescoping is only a slight extension of telescoping, in par- 
ticular, all the enhanced telescoping algorithms from Section [43] carry over to cre- 
ative telescoping. In all variations, a polynomial iTi;*-field (F, cj) (more precisely 
a depth-optimal HZ* -field 1671 for efficiency reasons) is constructed in which the 
summands F„+,(^) (0 < / < d) can be expressed. Starting from there, the following 
tactics are most useful to search for a solution of (1401 . They are activated by using 
the option Simplif yByExt->Mode where Mode is chosen as follows. 

• None: The solution G{k) is searched in (F,cj), i.e., only objects occurring in 
Fn+i{k) are used. Here a special instance of FPLDE is solved; see Example fTTI 

• MinDepth: The solution G{k) is searched in terms of sum extensions which are 
not more nested than the objects in F„^j{k) and which have minimal depth among 
all the possible choices ITMI . This is the default option. 

• DepthNumber: The solution is given in terms of sum extensions which are not 
more nested than Y!k=oPn+i{k), however, if the nesting depth is the same, the num- 
ber of the objects in the summands must be smaller than in F„+i{k). If such a recur- 
rence exists, the machinery from iJSZ) computes it. Using this refined version for our 
example, one finds a recurrence of order 1 (instead of 2) 

(;i + 1)^A_3(« + 1) + (« + 2)3A_3(n) 
= 6(« +2)(/i + \fSx{n) + {ln + U){n + \f + 3{n + 2fY.{n-2i) [\ (42) 



-3 



1=0 



where the sum E{n) ~ L!Lo(" ^ 2/)(") does not contain 5'i(/); it turns out that 
E{n) ~Q (using again our tools) and the recurrence simplifies further to 

(« + 1)^A_3(« + 1) + (« + 2)^A_3(«) = 6(« + 2)(« + l)^5i («) + (7« + 13)(« + 1)2. 

(43) 



4.3 Solving recurrence relations 

Next, we turn to recurrence solving in terms of nested hypergeometric sum expres- 
sions, i.e., expressions that evaluate to d' Alembertian sequence solutions. 

Example 12. Given the recurrence rec in Out|8] of A_3(«) = SUM [n] , all nested 
hypergeometric sum solutions are calculated with the following Sigma command; 

in[9]:= recSol = SolveRecurrence[rec,SUM[n]] 

Out[91. {{0.-(-l)°(ii+l)='},{0,(-l)-'(-Si(n)(n+l)=-(n+l)^)}. 

{l,6{n+ l)Si(n) + (-l)°{5(n+ l)=S_3(n) - 6(n+ Ij^S-s.iW) + 1}} 

The output means that we calculated two linearly independent solutions Hy (n) — 
-(-1)"(« + 1)^ and H2{n) = (-1)"( - 5i(«)(« + 1)^ - (« + l)^) (for « > 0) 
of the homogeneous version of the recurrence and a particular solution P(«) = 

6(n + l)5i(n) + (-l)"(5(« + l)35_3(n)-6(« + l)35_2,i(n)) + l (for « > 0) of 
the recurrence itself; since the solutions are indefinite nested, the verification of 
the correctness can be verified easily by rational function arithmetic. Note that 
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{ciHi{n) +C2H2{n) +P{n)\ci,C2 G Q} produces all sequence solutions whose en- 
tries are from Q. Since also A_j(n) is a solution of the recurrence, there is an el- 
ement in L that evaluates to A_3(«) for all n > 0. Using, e.g., the first two initial 
values A_3(0) = 1 and A_3(l) = 5 the ci,C2 are uniquely determined: ci ~ C2 = 0. 
Thus we arrive at A_3 («) = P{n), i.e., we discovered and proved the identity (|3]l for 
n > (recall that we verified that both sides satisfy the same recurrence and that 
both sides agree with the first two initial values). This last step is executed by taking 
recSol and mySum = A_3(«) (to get two initial values) as follows. 

In[i0]:= FindLinearCombination[recSol,mySum,n,2] 

Out[101= 6(n+ l)Si(n) + (-l)°(5(n+ l)^S_3(n) - 6(n+ l)^S_2.i(n)) + 1 

In general, Sigma can solve the following problem ISOl [121 l62l . 



Problem RS: Recurrence solving. 

Given polynomials ao{n),... ,ai„{n) G K{n) and a nested hypergeometric sum 
expression f{n). Find the full solution set of the mth-order linear recurrence 

ao{n)G{n) + ■ ■ ■ + a,„{n)G{n + m) = f{n) (44) 

in terms of nested hypergeometric sum expressions. I.e., return _L if there is no 
particular solution. Otherwise, find A G N and nested hypergeometric sum ex- 
pressions ((l,P(n)), (0,//i(n)),. .., (0, ///(«))) where P{n) is a particular solu- 
tion and Hi («), . . ■Hi{n) are solutions of the homogeneous version of (l44l i for 
« > A; the sequences (in ^(IK)) produced by H\ (n), . . ■Hi{n) are linearly inde- 
pendent. In addition, all sequences (G(«))„>o G K^, that are solutions of (l44li for 
all n > A and that can be given by nested hypergeometric sum expressions, can 
be also produced by 

L = {P{n)+ciHi{n) + . ..ciHi{n)\ci G K} (45) 

starting from « > A. 



The following solution rehes on ||39l l57l l62l |70l l76l . 
Solution*. Construct a polynomial J7Z*-field (F, a) as in Definition |23] with R de- 
fined as in ^ together with an evaluation function ev : /? x N — !■ K in which one 
obtains explicitly a G /? with A' G N such that ev(0,«) = /(n) for all n > A'; 
again Footnote |4] applies. In other words, with a,- := ai{x) G K(x) we can reformu- 
late (|44| with 

oo^ + aiGig) + ■■■ + a,„a"'{g) = (p. (46) 

Factorize the homogeneous recurrence (written as linear operator) as much as possi- 
ble in linear right factors using Hyper |l57| . Each linear factor describes a hypergeo- 
metric solution which is adjoined to our ni;*-field (see Problem RP); for simplicity 
we exclude the possible case that (—1)" is needed for this task. Applying Algo- 
rithm ll62l Alg. 4.5.3] to this recurrence returns the output (( 1, p), (0, /ii ),..., (0,/i/)) 
in a polynomial 7TZ*-field (E, a) that contains (F, a) with the following property: 
p is a particular solution of (l4&t and the /z, are / hnearly independent solutions of the 
homogeneous version of ( l46l l. We omit details here and remark only that it is crucial 
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to solve ( I22I 1 as subproblem. Then extend the evaluation function from F to E (we 
are free to choose appropriate lower bounds and constants of the sums/products), 
and let //,(«) (1 < / < /) and P{n) be the nested hypergeometric sum expressions 
that define the evaluations ev(/!,,«) and ev{p,n), respectively. Compute A such that 
the //,(«) are solutions of the homogeneous version and P{n) is a particular solution 
of (llll for all « > A; see Remark|3]with d = 0,fo := <P. 

Remark. (1) If one computes m = I linearly independent solutions plus a particular 
solution, the set ( |45] | gives all solutions. If this is not the case, the completeness of 
the method, i.e., that no solution in terms of nested hypergeometric sum expressions 
is missed, needs further justification: it can be deduced from Ii62i Cor. 4.5.2] and 
Remark[3j for deep insight and alternative proofs see ll36l and ll59l . 

(2) The derived solutions are highly nested: For each additional solution one needs 
one extra indefinite sum on top. In most examples the simplification of these solu- 
tions (see Subsection l4.1b is the most challenging task; see, e.g.. Example [l4l 

(3) Since the solutions are indefinite nested, the shifted versions can be expressed by 
the non-shifted versions. Using this property and considering the sums and products 
as variables, the correctness can be verified by rational function arithmetic. 

(4) Also the fl,(«) in ( |46] | can be from a i7Z*-field and one can factorize the differ- 
ence operator in linear factors; this is based on work by [E6] ITOl [TTl . 

Example 13. * We construct the polynomial iTr*-field {Q{x){h),a) with a{x) = 
x+1 and <y{h) = /1 + -^ over Q and interpret the elements with the evaluation func- 
tion ev : Q{x)[h] xN ^ Q{x)[h] canonically defined by ( |29] l and ew{h,n) = Si{n). 
In this way, we can reformulate the recurrence Out(8] with 

{x + 2)'{x + 3)^g + {x+l)\x + 3)H2x + 5)a{g) + {x + l)\x + 2fa\g)^ 

(20x^ + 138x2 + 311x + 229)(x+l)2 + 6(x + 2)2(x + 3)(2x + 5)(x+ 1)^/2. (47) 

Then we execute the recurrence solver in this 77Z* -field and get as output the differ- 
ence ring (Q(x)[m][/i][j'][//],c7) with a{m) = —m where m^ = 1, a{s) =s + 



and a{H) = H -\ iJ.W such that constaQ{x)[m][h][H] = Q. There it returns 



the linearly independent solutions hi = m{x+ ly and /z2 = m(^h{x+ !)■' + {x+ 1)^) 
of the homogeneous version of ( |47] | and the particular solution p ~ 6(ji:+ l)/j + 
m(5(x+ l)''^ -6(x+ 1)^//) + 1 of (021) itself. Note that the solutions (coming 
from the factorization of the recurrence) have been simplified already using the 
technologies presented in Section 14.11 Finally, we extend the evaluation func- 
tion from Q{x)[h] to Q{x)[m\[h][H] by ew{m,n) = (-1)", ev(s,«) = 5_3(n) and 
ev{H,n) = S-2.i(n). This choice yields H\{n) = ev(/ii,«), H2{n) = ev(/i2,n) and 
P{n) =ev{p,n) as given in Out[9]. 



Example 14. In ||23]| (see also B2ll ) recurrences are guessed with minimal order that 
contain as solutions the massless Wilson coefficients to 3-loop order for individual 
color coefficients 1481 . Afterwards the recurrences have been solved. The largest 
recurrence of order 35 could be factorized completely into linear factors in about 
1 day. This yields 35 linearly independent solutions in terms of sums up to nesting 
depth 34. Then their simplifications in terms of harmonic sums took 5 days. 
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5 Simplification of multiple sums with EvaluateMultiSums 

In Section |4] we transformed the definite sum A_3(«) to a nested hypergeometric 
sum expression given in (|3]l by calculating a recurrence and solving it. Applying 
this tactic iteratively leads to a successful method to transform certain classes of 
definite multiple sums to nested hypergeometric sum expressions. Consider, e.g., 

= :Fi(n.» 



n-2 ;+l (_-\\r(J+^\r\ -j+n+r-2 ( _-\\s (-j+n+r-2\ 

F{n):=j:{-j + n-2)lj: ) '\^^/] Y. \^\ -,/ (48) 

'^ . ' 

=:Fo{nJ,r) 

which arose in QCD calculations needed in f6\; see also fSl. We zoom into the sum 
Fo{n,j,r), a definite sum over a hypergeometric sequence. Calculating a recurrence, 
solving the recurrence, and combining the solutions leads to the simplification 

F„i„J,r)- ' . (-l)«0-+l)!(-;+"-l), 



(n + l)(-y + n + r-l) (n - l)n(n + l)(-y- l),(2-n)j 



This closed form could be also derived by hypergeometric summation 11581 : the dou- 
ble sum Fi («, j) turns out to be a single sum. Next, finding a recurrence for this sum 
and solving the recurrence lead to a nested hypergeometric sum expression w.rt. j: 

In other words, F (« ) can be written as a definite sum where the summand is a nested 
hypergeometric sum expression. Therefore we are again in the position to apply our 
technologies from Section|4] Computing a recurrence and solving it yields 

-«2_„_l (-l)"(„2+„ + l) Sijn) S2in) 2S^2{n) 
^"' «2(„ + i)3 + „2(„ + i)3 ^(„ + l)2 n + l n + l ' 

Summarizing, we transformed a definite nested sum from inside to outside to a 
nested hypergeometric sum expression. More generally, we deal with the following 



Problem EMS: EvaluateMultiSum. Given a definite multiple sum 

fim,n,k) 



L(){m.n)Li(m.n,k) Ly(m,n,k,ki,...,ky_i) 

F{m,n)^ Y^ ^ ... Y 7im,n,kM,...Ad (49) 

k=l ki=li ky=h. 

with a nested hypergeometric sum expression / w.rt. k^, integer parameters n and 
m = {m\ . . . ,mr), and L,(. . . ) being integer linear (see Footnote|5]l or °o. 

Find A G N and a nested hypergeometric sum expression F{m,n) such that 
F{m,n) = F{m,n) for« > A. 
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Method. Apply the techniques of Section|4]recursively as follows Il22l . 

1. Transform the outermost summand f{m,n,k) to a nested hypergeometric sum 
expression w.rt. k by applying the proposed method recursively to all the arising 
definite sums (i.e., the parameter vector m is replaced by (m, n) and the role of n 
is k). Note that the sums in / are simpler than F{m,n) (one definite sum less). If 
the summand / is free of sums, nothing has to be done. 

2. Solve Problem CT: Compute a recurrence (l3Fi for the sum A{n) = F{m,n); if 
this fails, ABORT. If successful (say it is of order o), the right hand side might 
be again an expression in terms of definite sums, but their summands are simpler 
than / (see, e.g., the recurrence (l42t). Apply the method recursively to these 
sums such that the right hand side is transformed to a nested hypergeometric 
sum expression w.rt. n (see, e.g., recurrence (|4Jt ). 

3. Solve Problem RS: Compute all nested hypergeometric sum solutions of the re- 
currence (l38T l and simplify the solutions using the techniques from Section |4T| 

4. Compute o initial values, i.e., specialize the parameter n to appropriate values 
from N, say n = 1,1 + l,...,l + o — I, and apply the method recursively to the 
arising sums where mi takes over the role of n and the remaining parameters are 
{1112,..., m,-). If no parameter is left, the expression is a constant. It is usually 
from Q (if no sum is left) or it simplifies, e.g., to multiple zeta values 1201 or 
infinite versions of S'-sums |9] and cyclotomic sums [HI . 

5. Try to combine the solutions to find a nested hypergeometric sum expression 
w.rt. « of F{in,n). If this fails, ABORT. Otherwise return the solution. 

Remark. (1) The existence of a recurrence in Step 2 is guaranteed in many cases (in 
particular for sums coming from Feynman integrals li24l ) by using arguments, e.g., 
form iHOl l58l |79l [TSl . Here often computation issues are a bottleneck. Usually, we 
succeed in finding recurrences when / consists of up to 100 nested hypergeometric 
sums. If / is more complicated (or if it seems appropriate), the sum is split into 
several parts and the method is applied to each sum separately. 

(2) Termination: The method is applied recursively to sums which are always sim- 
pler than the original sum (less summation quantifiers, less parameters, or less ob- 
jects in the summand). Hence eventually one arrives at the base case. 

(3) Success: If the method does not abort in one of the executions of step 2 or step 5, 
it terminates and outputs a nested hypergeometric sum expression w.rt. n. Note that 
finding not sufficiently many solutions of a given recurrence in step 5 is the main 
reason why the method might fail. For general multiple sums this failure would 
happen all over However, e.g, in the context of Feynman integrals, the recurrence 
is usually completely solvable (i.e., we find m linearly independent solutions of the 
homogeneous version of ( |38] | and one particular solution of the recurrence itself). 

We emphasize that 3-loop Feynman integrals with at most 1 mass 1241 can be 
transformed to multiple sums and that the simplification of these sums is covered 
exactly by Problem EMS. The described method is implemented in the following 
new package which uses the summation algorithms in Sigma; 

in[ii]:= << EvaluateMultiSums.m 

EvaluateMultiSums by Carsten Schneider © RISC 
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In addition it uses (some of the many) functions from J. Ablinger's package 
HarmonicSums lIZSl [78l [191 Isl [Tl l9l to transform -if possible- the arising indefi- 
nite sums to harmonic sums, 5-sums, cyclotomic sums or their infinite versions, to 
find algebraic relations among these sums, and to calculate asymptotic expansions 
of these sums for limit computations (this is needed if upper bounds in ( |49] l are °o). 

In[i2]:= << HarmonicSums. m 

HarmonicSums by Jakob Ablinger (c) RISC 

Then inserting the summand with the summation ranges of (l48T l and the information 
that there is the extra integer parameter n with 2 < « < oo we can activate the sim- 
plification of the sum ( |48] ) to a nested hypergeometric sum expression as follows. 

,n[i3,:. EvaluateMultiSum[ ' ,^^^/ ' > ' \„_,%,i/ ■ 



|-rii-ri;. ^11 — s;(s-ri.j 

{{s,0,-j + n + r-2},{r,0,j + l},{j,0,n-2}},{n},{2},{~}] 
-n+l) Si(n) S2(n) 2S_2(n) 
1)3 (ii+l)= n+l n+l 

Similarly, we can calculate the simplification given in identity (|4|: 

in[U]:=EvaluateMultiSum[(l-4(n-2k)Si(k))("j ,{{k,0.n}}.{n}.{0},{~}] 

„^,l,,j^ (10(n + l)Si(n)+3)(n+l) ^ (-l)°(^°)'^(n + 1)^ |^ 7 f (-IVCI) A (-l)^(^^)Si(i) 



-n=-n-l (-l)''(n=+n+l) Si(n) S2(n) 2S_2(n) 
°"""'° n2(n+l)3 '^ n2(n+l)3 "^(n+lj^^lT+T" n+l 



2n + 3 (4n(n+2)+3 

As mentioned already above, the multiple sums coming from many 2-loop and 3- 
loop Feynman integrals fit into the input class of the package EvaluateMulti- 
Sums. Here two extremes occurred: In ||6] |2lj about a million multiple sums 
(mostly triple and quadruple sums) were simplified. Using the package Sum- 
Production ||22 T| we merged the sums to several 100 basis sums where each of 
the summands required up to 20 MB memory. The other extreme are sums whose 
summands are in compact size, but the number of summations is large; one of the 
most complicated input sums from 11211 is, e.g., (|6]l. In both setups the transformed 
summands during the EvaluateMulti Sum method became rather large contain- 
ing complicated nested hypergeometric sums. Only in the last step these nasty sums 
vanished and the expected nice result popped up; note that already for the trans- 
formation of the sum ( |48] l this effect is visible. Summarizing, the summation al- 
gorithms based on enhanced difference field theory, presented in this article, were 
indispensable to master the challenging calculations as given, e.g., in ifTTl l6l l4l 1211 . 
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